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! ABSTRACT 

Aims. Precise photometric measurements of the upcoming space missions allow the size, mass, and density of satellites of 
exoplanets to be determined. Here we present such an analysis using the photometric transit timing variation (TTVp). 
Methods. We examined the light curve effects of both the transiting planet and its satellite. We define the photometric central 
, time of the transit that is equivalent to the transit of a fixed photocenter. This point orbits the barycenter, and leads to the 

■ photometric transit timing variations. 

' Results. The exact value of TTVp depends on the ratio of the density, the mass, and the size of the satellite and the planet. 
, Since two of those parameters are independent, a reliable estimation of the density ratio leads to an estimation of the size 
' and the mass of the exomoon. Upper estimations of the parameters are possible in the case when an upper limit of TTVp is 
known. In case the density ratio cannot be estimated reliably, we propose an approximation with assuming equal densities. 

■ The presented photocenter TTVp analysis predicts the size of the satellite better than the mass. We simulated transits 
' of the Earth-Moon system in front of the Sun. The estimated size and mass of the Moon are 0.020 Earth-mass and 0.274 

Earth-size if equal densities are assumed. This result is comparable to the real values within a factor of 2. If we include the 
r real density ratio (about 0.6), the results are 0.010 Earth-Mass and 0.253 Earth-size, which agree with the real values within 20%. 

' Key words. Planets and satellites: general - Methods: data analysis - Techniques: photometric 
' 1. Introduction 

' The discovery of planets in other solar systems ("exo- 
planets") led to the question of whether these planets 
also have satellites ("exomoons"). The photometric de- 
tection of such a satellite was first suggested by Sartoretti 
& Schneider (1999, SS99 in the foUowings; also discussed 
by Deeg, 2002), who gave a method for estimating the 
mass of this exomoon. The idea is to calculate the baricen- 
tric transit timing variations (TTVb). This quantity has 
been applied to the upper mass estimations of hypothet- 
ical satellites around transiting exoplanets (Charbonneau 
et al. 2005; Bakos et al. 2006; McCuUough et al. 2006; 
Steflten et al. 2005; Gillon et al. 2006). In this theory, the 
barycenter orbits the star with a constant velocity, and it 
transits strictly periodically. As the planet (and the satel- 
lite, too) revolves around the planet-satellite barycenter, 
the transit of the planet wobbles in time. The TTVb is 
based on finding the center between ingress/egress or from 
some other method known from eclipse minimum timing. 
The circular velocity around the star can be expressed via 
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the m* mass of the star and the Op semi-major axis of the 
planet as w = \J^^^- The maximal normal amplitude (i.e. 
half of the peak-to-peak amplitude) is finally 

TTVb = -B,, (1) 

V 

where Be = is the distance from the centerpoint 

of the planet to the barycenter. Here as is the semi- 
major axis of the satellite, and rris < rrip are the mass 
of the satellite and the planet, respectively. If TTVb is 
observed or if there is an upper limit, the mass is esti- 
mated by setting to the borders of Hill-stable regions, 

as — an = ap I j the Hill-radius. With this se- 
lection, rris remains the single parameter and it can be 
determined as 

w TTVb w TTVb , 2 \i/3 

ms,sS99 = TOp = [Snipm^) (2) 

as ap 

(Eq. (24) in SS99). This model can be refined by includ- 
ing the results of Barnes and O'Brien (2002), who sug- 
gests more stringent constraints on the survivability of 
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exomoons: the maximal distance from the planet is about 
an/S, depending slightly on the mass ratio. More recently, 
Domingos et al. (2006) have suggested that the maximum 
semimajor axis is 0.4895a// for a prograde satellite and is 
0.9309a/f for a retrograde satellite, both in circular orbit 
and with a mass ratio ms/rUp = 0.001. With the realis- 
tic choice of Os = aH/5, the satellite masses would lead 
to better estimates that are about three times larger than 
given by Eq. 2. However, this modification has been rarely 
adopted in the literature. 

In Szabo et al (2006) we designed numerical simula- 
tions to predict how probable the discovery of an exomoon 
is with a photometric technique. For this detection we sug- 
gested using the (photometric) central time of the transit, 
r as 



(3) 



where rrii is the i-th magnitude measurement at U time. 
We show later the photometric transit timing variation, 
TTVp calculated from this r requires its own physical 
interpretation, which differs from the conventional TTVb. 
The difference is to account for the photometric effects of 
both the planet and the satellite: even when the satellite 
cannot be detected directly in the light curve, its presence 
can actually dominate TTVp. We have shown that t is 
robust, therefore it is applicable even in the case of as 
"noisy" data series as the COROT and Kepler missions 
provide. These missions are able to detect a Moon-like 
satellite around an Earth-like planet with 20% probability 
(Szabo et al., 2006). 

In this paper we give the complete analytical descrip- 
tion of this theory of photometric timings, directly leading 
to estimation of various parameters such as the mass and 
the radius of the satellite. We test the formulae by mod- 
eling the Earth-Moon system. 

2. Light curve morphology 

Let us assume a planet and its satellite transiting the star. 
The position of the planet and the satellite is such that 
their images do not overlap as seen by the observer. We 
then assume that the relative position of the planet and 
the satellite does not change during the transit; the transit 
is central and the orbital inclination of the satellite is 0. 
Let the sampling of the light curve be uniform. Under 
these assumptions, r is exactly the barycenter of the light 
curve (i.e. considered as a polygon). Furthermore, let 



fi^ms/nip, x^Ps/Pp, ^^rs/rp 



(4) 



be the ratio of the masses, densities and sizes (radii), re- 
spectively. Any of the three can be expressed from the 
two other ones. In the case of the direct detection of the 
satellite in the shape of the light curve, one can directly 
calculate = Arris/Amp, where Am^ and Arup are the 
brightness decreases due to the satellite and the planet, 
respectively. But even if the noise level is too high for a 
direct detection, the presence of a moon can be reasoned 



Am, 



Satellite alone 
+ 

Planet alone / 




Together 



Fig. 1. Composition of the transit light curve from Amp 
and Am = . 



indirectly (Szabo et al. 2006), as formulated in the follow- 
ing. 

The shape of the light curve is the sum of two compo- 
nents: a single planet and a single moon in transit (Fig. 
1). The shape of those components is 



Amp^ti = — f{ti - Tp), 



(5) 



(6) 



where Amg^u and Arup^u are the magnitude decreases at 
time ti, and Tp and Ts are the times when the planet and 
the satellite each passes alone before the central meridian 
of the star. Here f{x) is the normalized shape function 
of the transit. We can assume that f{x) is an axially 
symmetric light curve to a; = 0, and its off-transit value 
is 0. 

If we consider the transit of the planet and the moon 
together. 



Ami — AvTLs^ti + Amp^t;. 



In this case we know from Eqs. 3 and 7 that 



Amj 



(7) 



(8) 



Let N — ^f{ti) be the "area" under the function /. 
With this, we can easily rewrite the nominator as 



(9) 



Taking the axial symmetry of / into account, the numer- 
ator of Eq. 8 can be similarly rewritten as 



t.ArUsM + V t.Amp^ti = -^N ■ r, + -|iV • Tp, (10) 
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and finally 



T = 



- TpTp 



(11) 



Now let us consider an extreme configiiration with a 
leading moon. We consider the first contact time, when 
the ingress of the satellite begins, and relate it to the time 
when the satellite is exactly in the center of the stellar 
disc. Between these events, Ts = -^{rs + J^*) time has to 
elapse. The time between the first contact of the satellite 
and the time when the planet is in front of the center of 
the stellat disk is Tp = ^{Vg + r^, + Us). Substituting these 
and the analogous formulae into Eq. 11, 



V \ ri+ri. 



+ r* 



with a leading satellite and 



+ r* 



rp2 _|_ rp2, 



(12) 



(13) 



with a trailing satellite. 

The basic conclusion is that the planct-satcUitc sys- 
tem can be substituted by a single celestial body, which 



can cause the same TTVp and 
planet-moon line at Pc = 



which is located on the 

„2 

distance from the 



r -\-r 



planet. We call this particular point as "photocenter" in 
the following. 

The barycenter lies at Be = „°°"^ distance from the 
planet and it revolves strictly periodically. The photocen- 
ter lies somewhere else between the planet and the satel- 
lite, therefore the r time of the transit will wobble in time 
(See Fig. 2). The time between the transit of the barycen- 
ter and the photocenter is TTVp, which can be expressed 
as 



TTV, 



\Pa - BJ = 



(14) 



With neglecting -d"^ and /i to 1 (i?^ ^ 1, <C 1, their 
values ^ 0.01) in the nominators, TTVp simplifies to our 
basic equation 



TTVp = ^ 

V ' 



(15) 



where /U = mp/nis is the mass ratio and •& = Vs/vp the 
ratio of radii. The comparison of Eqs. 1 and 14 shows the 
difference between the SS99-thcory and the present ap- 
proach: the introduction of the Pc term. This reflects the 
different basic concepts. In the SS99-theory, TTVb is a 
dynamical timing effect caused by the revolution of the 
planet around the barycenter. In our present approach, 
TTVp is due to the revolution of the photocenter, which 
both combines the dynamical and the photomctrical prop- 
erties of the moon and offers a potential to estimate the 
mass and the size of the satellite. 




as-Pc+r, 




Fig. 2. Position of the barycenter and the photocenter in 
a transiting system (not to scale). 
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Fig. 3. The dependence of TTVp (multiplied by the term 
v/cs) on the size ratio. The three different curves show 
three satellites with different density ratios. 
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Fig. 4. The same as Fig. 3 but for the mass dependence. 
Note that TTVp depends noticeably on the density in this 
model, which makes the mass determination more ambigu- 
ous. 
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2.1. Expressions for the radius and mass of tine moon 

By expressing the mass ratio /i through the radii r^, 
and the densities pp, ps in Eq. [151 we get the formula for 
the size ratio d: 



TTVp - 



Ps_ 

Pp 



V 



ATTVp = — { [2i9 - 3d^x\ Ai9 - ^^Ax} 



(16) 



(17) 



where x ^ Ps/Pp is the ratio of the densities. Then, 
ATTVp expresses the error propagation from -d and x 
to TTVp using the total derivative of Eq. 16. We note 
that in most cases — X > and the absolute value can 
be replaced by normal brackets. This is because (i) in the 
case of a giant planet << 1, x > 1 can be expected, and 
(ii) in the case of an Earth-like planet d < 1 and x ~ 1 is 
plausible. If is known from a direct photometric detec- 
tion, this equation can be used to determine the density 
ratio. 

From Eq. [T5] we can eliminate the radii using masses 
and densities. This gives the formula for determining the 
mass ratio. 



TTVp = — 

V 



ATTV:, 



nip/pp 



2 








3 


rris 


as 














mp 


V 





All ~ ^Ax 

3x3 



(18) 



(19) 



The error propagation analysis shows that Eq. [16] is quite 
stable against the ambiguity in X: while Eq. [18] is more 
sensitive to x- This means that the photometric method is 
better for a size determination than for a mass determi- 
nation. This is also demonstrated in Figs. 3 and 4, which 
plot TTVp using Eqs. [16] and [18] for different sizes and 
masses for some discrete values of x- 

2.2. X = 1 estimations 

In the case of the known transits, the x relative density 
is unknown. To make an order-of-magnitude estimate or 
an upper limit of the size and the mass of a hypothetical 
satellite, we have to choose some value for the density 
ratio, e.g. x = 1- This choice leads to 



TTVn 



-j?2(i_^) 



TTVp « — 

V 



2/3 



(20) 



(21) 



If we assume that = a/z/S, these equations contain 
only one unknown parameter. As the size determination 
is not too density-dependent (Eqs. 17, 19., Figs. 3, 4), one 
can expect better results using Eq. 20 for size determina- 
tion, but Eq. 21 can also be used for a very rough mass 
estimation. 
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Fig. 5. Grid representation of TTVp caused by a satellite 
with different x density ratios and size ratios. 



2.3. Testing with the Earth-Moon system 

A simple test of the calculations described above is to ex- 
amine the Earth-Moon system. In Szabo et al. (2006) the 
corresponding numerical calculations are described with 
the real sizes, orbital periods, and the solar limb darkening 
also taken into account. The resulting TTVp is 15.2 min- 
utes (normal amplitude). The corresponding TTVpu/ag is 
0.054, when substituting as = a_ff/3. 

First we calculated a x = 1 models leading to the same 
timing effect. We found that ■& = 0.274, ^ = 0.020. The 
size of this Moon-model agrees well with the real size of 
the Moon, and the mass is about twice the real value. This 
precision is, however, acceptable for a trial estimation for 
the parameters of an exomoon. However, if we somehow 
got to know that x = 0.605 in the Earth-Moon system, we 
would get an acceptably precise result for both parameters 
asd = 0.253, p. = 0.010 (Table 1, 2nd row). Both values 
are concordant with the real parameters within 20%, and 
suggest the reliability of this method. We conclude that 
he determined sizes agree with the real properties, and 
an order-of-magnitude estimate for the mass of the Moon 
is also possible even when no information on the density 
assumed. 

In order to compare the results x = 1 estimations to 
other possible density ratios, a grid presentation of the 
timing effect is plotted in Fig. 5 based on Eqs. 15 and 
18. The lines of constant x and of constant d values are 
plotted in the ^ vs. TTVp -v/as space. With this grid one 
can determine how the mass and the size estimates of the 
Moon vary if we have a reliable estimation for x- Also, 
this grid may help in estimating the size and the mass 
of a possible exomoon in the case of the positive TTVp 
detections of the future. 
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l/3aH TTVp TTVp.^;/a, ^ = = ^ 
(km/s) (10^ km) (s) 

X = 1 29^8 499^6 912 01)54 (L274 01)20 

X = 0.605 29.8 499.6 912 0.054 0.253 0.010 

real 0.272 0.012 

Table 1. Test analysis of the simulated transit light curve of the Earth- Moon system, using data sets with the noise 
and sampling properties of the Kepler mission. 



3. Conclusions 

We have presented a novel approach to timing-effect mod- 
eling of an cxoplanct with a satellite. This approach of- 
fers easy calculations with the light curves and helps es- 
timate the mass and the size of a possible satellite. The 
application to high-precision photometric data, such as is 
provided by the upcoming satellite missions will lead to 
estimating the masses and sizes of the moons. The differ- 
ence between this approach and the SS99-theory is that 
here we have defined the photometric central time of the 
transit r, which accounts for the light curve distortions 
due to the planet and the satellite, too. Although these 
seem to be tiny corrections, they have a significant effect 
in the interpretation of the measurements. An important 
consequence is that TTVp cannot exceed a limit if the 
mass (or the size) of the satellite is increased more and 
more. To illustrate this, let us imagine a "double planet", 
a planet and a satellite having exactly the same size and 
mass. Because the positions of this "planet" and "satel- 
lite" are persistently symmetric, one suspects no timing 
effect to occur. Our formulation reproduces this result, if 
we include i? = /z = 1 in Eq. 15. 

In Fig. 5, the grid representation of the TTVp leads 
to similar conclusions. The constant % lines at low fj, val- 
ues shows that the bigger the size of a moon, the larger 
the TTVp result. But this effect does not increase con- 
stantly with size, e.g. a satellite with x = 3 displays maxi- 
mal TTVp if = 0.034, and then TTVp starts decreasing 
with increasing ji. The constant i9 lines, on the other hand, 
show that TTVp due to a satellite of a fixed size is maxi- 
mal when the density is low. This is because the barycen- 
ter in this system gets closer to the center of the planet, 
while the photocenter remains at its position regardless 
of X- Therefore the distance between the photocenter and 
the barycenter increases, thus TTVp increases as well, but 
never exceeds the time between the transit of the photo- 
center and the center of the planet. 

The limiting values of TTVp referring to low-density 
satellites can be calculated from Eqs. 16 and 18 gener- 
ally and from Eqs. 20 and 21 in the x = 1 case, and 
the highest possible value of TTVp can be calculated. In 
this case the timing effect is the largest with the value 
of wTTVp/tts w 0.296 if = 2/3 or, what is the same, 
/i = 8/27. From the viewpoint of the observation, too large 
an upper estimation of TTVp is meaningless whenever it 
exceeds the above limits. 



However, a higher value of TTVp may also be inter- 
preted physically. There are more processes that result in 
timing effects that can exceed the TTVp caused by an 
exomoon. A well-known example is the perturbation of a 
second planet (Agol & Steffen 2005; Steffen 2006; Gillon 
et al. 2006; Heyl & Gladman 2006). Ford & Gaudi (2006) 
suggest that the presence of massive exotrojan asteroids 
can lead to the shift in the transit central time regard- 
ing the time of zero radial velocity difference between the 
planet and the barycenter. The libration, even on a horse- 
shoe orbit, may result in a timing effect that can exceed 
a few minutes. 

The results of this paper can be summarized as follows. 

— We gave a theoretical framework for a photometric 
timing effect caused by moons of extrasolar planets 
during transits. We included the photometric processes 
due to the moon via summing the entire photometric 
signal. This concept offers more parameters to be de- 
termined, such as the sizes and/or the masses of the 
satellites. 

— The results are supported by the numerical simulations 
of Szabo et al., and another simple test is described 
here. With the analysis of the simulated transits of 
the Earth-Moon system in front of the Sun, we could 
almost reproduce the parameters of our Moon exactly. 

— We argued that the size of an exomoon is a better 
predictable parameter than its mass, and we suggest 
using this estimation in further analyses. 
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